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While supernova remnants (SNRs) have long been considered prime candidates for the source 
of cosmic rays, at least to energies up to ~ 10 14 eV, it is only over the past several years that 
direct evidence of such energetic particles in SNRs has been uncovered. X-ray observations of 
several shell-type SNRs have now revealed sites dominated by nonthermal emission, indicating 
an electron population whose energy extends far beyond the thermal distribution typical of 
such SNRs. In other remnants, discrepancies between the shock velocity and the electron 
temperature points to a strong cosmic ray component that has essentially thrived at the 
expense of the thermal component of the gas. Modeling of the radio, X-ray, and gamma-ray 
emission provides strong constraints on the acceleration mechanism as well as the properties 
of the ambient medium in which the mechanism prospers. In the innermost regions of some 
SNRs, particle acceleration is taking place over much different scales. The formation of Crab- 
like pulsar wind nebulae (PWNe) is understood to require the presence of a termination shock 
at which the relativistic pulsar wind is forced to join the slow expansion of the outer nebula. 
While the acceleration mechanism is necessarily different, these shocks also act as sites in which 
particles are boosted to high energies. In the Crab Nebula, optical wisps mark the location of 
this termination shock. Recent X-ray observations have begun to reveal the termination shock 
zones in other PWNe, and are now allowing us to constrain the nature of the pulsar wind as 
well as the flow conditions in the outer nebula. Here I present a summary of the properties 
of shock acceleration in these two distinct regions of SNRs, and review recent observational 
results in which the properties of the shocks are finally being revealed. 



1 Introduction 



Supernova remnants and their associated pulsars are sites in which strong shocks act to accelerate 
particles to extremely high energies. The connection between SNRs and the energetic cosmic 
rays that pervade the Galaxy has long been assumed, for example; shock acceleration by the 
SNR blast wave provides ample energy for the production of multi-TeV particles, and recent 
observations of nonthermal X-ray and VHE 7-ray emission from several SNRs has confirmed the 
presence of electrons at these high energies. At the same time, models for the structure of PWNe 
predict particle acceleration at the wind termination shock, and recent X-ray observations have 
begun to probe these acceleration sites. Here I review recent and ongoing observational X-ray 
studies of particle acceleration by shocks in SNRs and PWNe. 



2 Shocks and Particle Acceleration in SNRs 

The birthrate and overall energetics of SNRs provide a strong plausibility argument that they 
are a major source of cosmic rays up to the "knee" of the spectrum. However, until recently 



the only direct evidence of energetic particles from SNRs has been from the radio synchrotron 
emission. But this corresponds to electron energies 
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where v is the frequency of the radio emission and is the magnetic field strength in micro- 
Gauss. Thus, while providing a clue, the radio emission samples only electrons far below the 
cosmic ray knee. X-ray observations allow us to characterize the thermodynamic states of SNRs 
and to infer the dynamics of their evolution, probing much higher particle energies. As the blast 
wave from a supernova explosion expands, material from the surrounding circumstellar material 
(CSM) and ISM is swept up into a shell of hot gas. For an ideal gas, the shock jump conditions 
yield an increase in density by a factor of 4 and a postshock temperature of 
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where V s is the shock speed, m is the proton mass, and \i is the mean molecular weight of 
the gas {fx ~ 0.6). This shock-heated gas yields the familiar X-ray emission, characterized by a 
thermal bremsstrahlung spectrum accompanied by strong emission lines. As increasing amounts 
of material are swept up by the blast wave, the shock is decelerated and the ejecta from the 
explosion encounter the dense shell. A reverse-shock is generated, heating the ejecta. At early 
times, then, the X-ray spectrum is dominated by emission from the ejecta. As the amount of 
swept-up material increases, the spectrum becomes dominated by emission from material with 
ISM abundances. 

In addition to thermal heating of the swept-up gas, some fraction of the shock energy density 
goes into nonthermal production of relativistic particles through diffusive shock acceleration. 
The maximum particle energy in such a scenario can be limited by radiative losses, the finite 
age of the SNR, or particle escape from the accelerating region. Radio observations of SNRs 
provide ample evidence of electrons with GeV energies through synchrotron radiation in the 
compressed magnetic field at the remnant shell. At higher particle energies, 7-ray production 
may result from nonthermal bremsstrahlung of electrons colliding with ambient gas, from inverse 
Compton scattering of ambient photons, and from the decay of neutral pions created by the 
collision of energetic protons. If the relativistic particle component of the energy density becomes 
comparable to that of the thermal component, the shock acceleration process can become highly 
nonlinear. The gas becomes more compressible, which results in a higher density and enhanced 
acceleration. When the acceleration is very efficient, the relationship between the shock velocity 
and the mean postshock temperature is no longer well approximated by Eq. 1; the acceleration 
process depletes thermal energy, and the temperature for a given shock velocity drops below 
that expected in the test particle casea. This process has been considered in detail by Baring 
et al. 13 who present a model for the radio to 7-ray emission. The broadband spectra of SNRs 
depend highly on ambient conditions, and X-ray studies of these SNRs reveal these conditions 
and can provide spectral measurements which strongly constrain the models. 

In the simplest picture, the passage of material through the SNR shock results in electrons 
and ions being boosted to the velocity of the shock. Because of the mass difference, this means 
that the electrons and ions are not initially in temperature equilibrium. The maximum timescale 
for equilibration is that provided by Coulomb interactions, but plasma processes may reduce 
this considerably. The state of equilibration is important, because while the dynamics of SNR 
evolution are dominated by the ions (which carry the bulk of the momentum) , it is the electrons 
that produce the X-ray emission we observe. Thus, when temperature measurements are used 
to infer the shock velocity, for example, the state of temperature equilibration is exceedingly 
important. Chandra HETG observations of SN 1987A,H for example, yield spectra consistent 



with an electron temperature of 2.6 keV. However, the broadened line profiles indicate a blast 
wave speed of ~ 3500 km s _1 which corresponds to a post-shock temperature of 17 keV, pro- 
viding evidence for incomplete electron-ion temperature equilibration. Similarly, the blast wave 
speed inferred from an expansion study comparing a Chandra image of IE 0102.2-7219 with high 
resolution images taken with Einstein and R OSAT indicates a post-shock temperature which is 
much higher than the observed electron temperature.! In this case, the discrepancy appears to be 
larger than can be accounted for assuming Coulomb equilibration, suggesting that a significant 
fraction of the shock energy has gone into cosmic ray acceleration rather than thermal heating 
of the postshock gas. Higher resolution X-ray expansion studies of IE 0102.2-7219 (and other 
young SNRs) are needed to confirm this scenario, but the notion that particle acceleration is 
efficient enough in some SNRs to significantly affect their dynamics has strong foundations from 
recent studies of other SNRs. As we discuss below, direct evidence of very energetic electrons 
now exists for a handful of shell-type SNRs. In addition to the three SNRs which we discuss in 
the following sections, for which synchrotron radiation dominates the X-ray emission, evidence 
for energetic particlesin the form of hard tails in the X-ray spectra have also been observed for 
Cas AH RCW 86,11 and other SNRs. 

2.1 SN 1006 

The first evidence of multi-TeV particles directly associated with a shell-type SNR was uncov- 
ered in studies of SN 1006 with the ASCA observatory! While spectra from the central regions 
of the SNR show distinct line emission associated with shock-heated gas, the emission from the 
bright limbs of the remnant was shown to be completely dominated by synchrotron emission. 
Reynolds E-3 modeled the emission as the result of diffusive shock acceleration in a low den- 
sity medium with the magnetic field orientation providing the distinct "bilateral" morphology. 
Subsequent observations with the CANGAROO telescope revealed VHE 7-ray emission from 
one limb of the SNR, Ej thus confirming the presence of extremely energetic particles. The low 
ambient density of this remnant, which resides well above the Galactic Plane, is insufficient 
to explain the 7-ray emission as the result of ir° decay from proton-proton collisions. Rather, 
the emission results from inverse-Compton scattering of microwave background photons off the 
energetic electron population in SN 1006. Using joint spectral fits to the radio as well as thermal 
and nonthermal X-ray emission, Dyer et al.@ conclude that the total energy in relativistic par- 
ticles is ~ 100 times the energy in the magnetic field, confirming the notion that the nonthermal 
particle component contributes significantly to the dynamics of the blast wave evolution. 

2.2 G 3 47. 3-0. 5 (RX J 17 13. 7-3946) 

ASCA observations of G347.3-0. 5§§ established this as the second member of the class of shell- 
type SNRs for which the X-ray flux is dominated by synchrotron radiation. Unlike SN 1006, this 
SNR appears to have evolved in the vicinity of dense molecular clouds with which the shock 
may now be interacting. CANGAROO observations E-3 reveal VHE 7-ray emission from the 
vicinity of the northwest rim, which is brightest in X-rays. Combining radio measurements from 
the ATCA with the X-ray and 7-ray results, Ellison, Slane, & Gaenslerlll used diffusive shock 
acceleration models to conclude that the radio and X-ray emission results from synchrotron 
radiation from a nonthermal electron population accelerated by the forward shock, and that 
the 7-ray emission can be self-consistently modeled as inverse-Compton emission (Figure 1). 
Combined with limits on the ambient density based on the lack of thermal X-ray emission, the 
models indicate very efficient particle acceleration with > 25% of the shock kinetic energy going 
into relativistic ions. A comparison of Chandra observations of the northwest rim with high 
resolution radio maps from the ATCA shows good overall agreement with the radio and X-ray 
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Figure 1: Comparison of broadband emission from the northwest limb of G347.3-0.5 with diffusive shock accel- 
eration model results (from Ellison, Slane, & Gaensler 2001). The radio and X-ray emission is from synchrotron 
radiation while the VHE 7-ray emission is from inverse-Compton scattering. Sensitivities of future high energy 

telescopes are shown for comparison. 



morphology in this region (Figure 2Y consistent with the interpretation that the emission comes 
from the same electron population. £3 

The nearby unidentified EGRET source 3EG J1714-3857 has been suggested as being as- 
sociated with G347.3-0.5, possibly resulting from the decay ofjieutral pions produced in the 
collision of accelerated ions with the nearby molecular clouds, Ii3 although this would appear 
inconsistent with the X-ray measurements unless the emission originates from a distinct spatial 
region. Most recently, Enomoto et al. present new CANGAROO data which, they argue, 
establishes 7r°-decay as the mechanism by which the TeV 7-rays are produced. However, the 
predicted spectrum overpredicts emission from the EGRET band by a large margin, @ making 
the claim for direct observation of ion acceleration appear problematic. 




Figure 2: Teft: ROSAT PSPC image of G347.3-0.5. Right: Chandra image of northwest rim of G347.3-0.5, 

with radio contours from ATCA. 



2. 3 G266. 2-1 . 2 (RX J 0852. 0-4 622 ) 

G266.2-1.2 is an SNR lying along the line of sight to the southeast corner of the Vela SNR. 
Discovered in the ROSAT All-Sky SurveyEI the remnant is evident only at higher X-ray energies 
where the very soft spectrum of Vela does not contribute. ASCA observations!^ show that the 



X-ray emission from this remnant is also predominantly synchrotron radiation. To date, no high 
energy 7-ray emission has been reported for this remnant, although there have been claims of 
44 Ti line emissions which, given the very short half-life, would require an extremely young age. 
The radio flux from G266.2-1.2 is extremely low, El much like that of SN 1006 and G347.3-0.5. 
The overall properties of this remnant are indicative of efficient shock acceleration of particles 
to very high energies, but further investigations are required to quantify the dynamics. 

3 Energetic Particles from Pulsars and Their Wind Nebulae 

Our basic understanding of "Crab- like" SNRs stems, iom the picture presented by Rees & 
Gunned , and expanded upon by Kennel & Coroniti,c3c3 in which an energetic wind is injected 
from a pulsar into its surroundings. As illustrated schematically in Figure 3, the basic structure 
of a pulsar wind nebula is regulated by the input power from the pulsar and the density of 
the medium into which the nebula expands; the pulsar wind inflates a magnetic bubble which is 
confined in the outer regions by the expanding shell of ejecta or interstellar material swept up by 
the SNR blast wave. The boundary condition established by the expansion at the nebula radius 
r^r results in the formation of a wind termination shock at which the highly relativistic pulsar 
wind is decelerated to v ~ c/3 in the postshock region, ultimately merging with the particle 
flow in the nebula. The shock forms at the radius r w at which the ram pressure of the wind is 
balanced by the internal pressure of the pulsar wind nebula: 

r 2 w = E/^rjcp) (3) 

where E is the rate at which the pulsar injects energy into the wind, r/ is the fraction of a 
spherical surface covered by the wind, and p is the total pressure outside the shock. 




Figure 3: Schematic view of a pulsar and its wind nebula. See the text for a complete description. (Note the 
logarithmic size scaling in the PWN figure when comparing with images shown elsewhere in the text.) 

The dynamics of the particle flow yield 7 ~ 10 6 for the electrons in the postshock region. H 
However, for typical magnetic field strengths the observed X-ray emission requires 7 > 10 8 . 
Particle acceleration at the termination shock apparently boosts the energies of the wind particles 
by a factor of 100 or more, to energies in excess of ~ 50 TeV. Arons & Tavani c3 (see also 
AronsH) note that this process cannot proceed by normal diffusive shock acceleration because 
the magnetic field at the termination shock must be nearly perpendicular to the flow. Rather, 
they argue that the acceleration is the result of resonant cyclotron absorption of low frequency 
electromagnetic waves emitted by ions gyrating in the compressed -B-field of the hot post-shock 
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Figure 4: Chandra image of the Crab Nebula (left) and the pulsar wind nebula associated with PSR B1509-58 

(right; image courtesy of B. Gaensler). 

gas. Ultimately, the pressure in the nebula is believed to reach the equipartition value; a 
reasonable pressure estimate can be obtained by integrating the radio spectrum of the nebula, 
using standard synchrotron emission expressions, and assuming equipartition between particles 
and the magnetic field. Typical values yield termination shock radii of order 0.1 pc, which yields 
an angular size of several arcsec at distances of a few kpc. 

The magnetic field in the PWN builds up with radius as the result of particle currents and 
wound-up magnetic flux from beyond the pulsar light cylinder at rj_,c = c/Sl. The finite lifetime 
of the synchrotron-emitting particles in the nebula, 

t syn » 5 x uP/{EV*vW) s 

(where v is the frequency of the synchrotron radiation), results in fewer high energy particles 
at large radii. This manifests itself as a variation of the nebula radius with energy (which is 
clearly seen in the Crab Nebula, for which the radio size is larger than that seen in X-rays) or, 
equivalently, as a radial variation in the X-ray spectrum, with the power law index steepening 
with radius. 

The pulsar wind can be characterized in terms of E and the parameter parameter a repre- 
senting the ratio of Poynting flux to particle flux. For the Crab Nebula, Kennel & Coroniti H 
find that small values of a (~ 10 -3 ) are required, indicating a particle-dominated wind. Yet 
current understanding of pulsar outflows predicts a 1 where the wind is launched, i! Some- 
how, between the pulsar light cylinder and the wind termination shock the balance between the 
electromagnetic energy and the kinetic energy of the flow changes dramatically. The ability to 
identify the termination shock and measure the emission parameters in the immediate region 
is now providing constraints on the nature of the wind that may ultimately help unravel this 
problem. Below I summarize recent X-ray investigations that have finally begun to probe this 
important shock region in which particle acceleration is apparently taking place. 

3.1 Crab Nebula 

The Crab Nebula is the best known of the class of pulsar wind nebulae and has inspired much 
of the theoretical work on PWNe. Powered by an energetic central pulsar, it emits synchrotron 
radiation from radio wavelengths up beyond the hard X-ray band. Optical wisps are observed 
in the inner nebula, at a position interpreted as the pulsar wind termination shock, EH and high 
resolution X-ray observations (Figure 4, left) reveal a distinct ring of emission in this same 
region as well as a jet emanating from the pulsar. Moreover, monitoring observations of the 



I -10"3«>" 

i 

-3i'30" 



la^ss' ai s as 1 32= 10 20 30 40 

Radius (arcsec) 

Right Ascension (JSOOO) 

Figure 5: Left: Chandra image of G21.5-0.9. The inset shows the brightness profile of the compact central region 
compared with a histogram of that from a point source in the telescope. Right: Spectral index variation with 
radius in the nebula, indicating the synchrotron burn-off of particles injected at the center of the nebula. 

Nebula 0E1 show that these and other detailed features are highly dynamic. The discovery of 
radio wisps in inner ring region 11 suggests that the acceleration site may be the same for the 
entire population of electrons that produce the broad-band synchrotron emission. 

In addition to the jet and inner ring, the X-ray image reveals an outer toroidal structure 
that presumably lies in the equatorial plane, as well as fine structure correlated with optical 
polarization measuremente, indicating that the structures trace the magnetic field. The early 
models of Rees &i Gunn E3 and Kennel & Coroniti I 26 !! 27 ! predict these basic properties as the 
result of a wound-up magnetic field, the large-scale confinement of the Nebula by the (unseen) 
supernova ejecta, and the termination of the pulsar wind flow by an inner shock. This picture 
leads to the inference of a low-cr wind described above. As we describe below, recent observations 
of PWNe have begun to show many similar features, indicating that our basic picture - while 
still poorly understood in detail - can at least be said to apply to a "class" of objects. 

3.2 PSR B 1509- 58 

The PWN powered by the energetic young pulsar B1509-58 (Figure 4. right) displays a complex 
morphology rivaling that of the Crab Nebula. Chandra observations Ell reveal a clear symmetry 
axis, presumably corresponding to the pulsar spin axis, with an apparent curved jet along this 
axis. Arc-like features in the inner nebula appear to correspond to equatorial structures similar 
in nature and origin to the wisps seen inside the Crab X-ray torus. Gaensler et al.B show that 
a ~ 0.005 at these wisps, with other compact knots seen very close to the pulsar indicating 
even smaller a values upstream of the wind shock zone. The PWN has a much weaker magnetic 
field than the Crab, and synchrotron losses affect the spectra of outflowing particles only at 
large radii. The hard spectra observed for the compact knots suggest a different acceleration 
mechanism than for larger features, downstream of the termination shock. In the northernmost 
regions of the PWN, clumps of thermally-emitting gas from the associated SNR appear to be 
embedded in the nebula. 

3.3 G 2 1.5-0. 9 

G21.5-0.9 is a compact SNR that exhibits strong linear polarization, a flat spectrum, and 
centrally peaked emission in the radio band. The SNR has a lower L x /L r ratio than the Crab; 
it is a factor of ~ 9 less luminous in the radio and a factor of ~ 100 less in X-rays. To date, 
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there has been no detection of a central pulsar. However, Chandra observations by Slane et 
al. H reveal a compact source of emission at the center of the remnant as well as a radial 
steepening of the spectral index consistent with synchrotron burn-off of high energy electrons 
injected from a central source (Figure 5). Using an empirical relationship between the total X-ray 
luminosity of the PWN with the spin-down power of the pulsar powering the nebula suggests 
the presence of a pulsar with E = 10 37 ' 5 erg s _1 , although Chevalier argues that the spectral 
variations in the nebula imply more efficient conversion of E into X-ray emission, and suggests 
that E ~ 10 36 ' 7 erg s _1 is more likely for the pulsar in this PWN. Detection of pulsations from 
the central source are required to address this further. Using the larger E estimate along with 
pressure estimates from the radio spectrum, Eq. 3 predicts a wind termination shock at a radius 
of ~ 1.5ry -1 / 2 arcsec from the pulsar, assuming a distance of 5 kpc. As indicated in the inset 
to Figure 5, the brightness profile of the compact X-ray source in G21.5-0.9 is broader than 
that for a point source. The ~ 2" extent of the source is consistent with the expected size of 
the termination shock zone. The small radius of the termination shock combined with the large 
radius of the PWN yields a <C 1; as with the Crab, the pulsar wind is particle-dominated in the 
termination shock zone. 

Slane et al. H also report the presence of a faint extended shell of emission surrounding 
G21.5-0.9 whose featureless spectrum may be associated with energetic particles accelerated 



by the SNR blast wave, similar to those observe 
Section 2. Safi-Harb et al.E3 and Warwick et al.t 



for SN 1006 and other SNRs as described in 
confirm the presence of the shell, but suggest 



that it is a faint extension of the PWN, though no extended radio shell is observed. 



3.4 3C58 

3C58 (Figure 6, left) is a young Crab-like supernova remnant. Historical evidence strongly 
suggests an association of the remnant with supernova SN 1181, which would make 3C58 younger 
than the Crab Nebula. Recent Chandra observations have identified the young 65 ms pulsar 
J0205+6449 at its center, e3 embedded in a compact nebula which appears to be confined by 
the pulsar wind termination shock. Ell The central region of 3C58 is shown in Figure 5 (center). 
The emission is clearly extended, with elongation in the N-S direction, perpendicular to the long 
axis of the main nebula; the full extent in this direction is ~ 25". PSR J0205+6449 resides at 
the center of this compact nebula. 

The western edge of the compact nebula lies directly along a radio filament, El shown as con- 
tours in Figure 5 (center!, and there appears to be a slight flattening of compact nebula along 
this side. Frail & MoffettEj suggested that this filament may represent the position of the termi- 
nation shock where the pulsar wind is confined by the interior pressure of the PWN. Adopting 
the generally accepted distance d = 3.2c?3.2 kpc, integration of the radio synchrotron spectrum 
yields a pressure of P ne b = 3.2 x 10 _10 <i£ 2 dyne cm~ 2 under the assumption of equipartition 
between the electron and magnetic energy densities. The spin-down properties of the pulsar 
give E = 2.6 x 10 37 /45 er g s ~ 1 where I45 is the NS moment of inertia in units of 10 45 gm cm 2 . 
Ram pressure balance should thus occur at r w = 5.5 x lO 17 /^ 2 r}~ x l 2 d\ % cm, or at an angular 
distance 9 = II.4I45 77"" 1 / 2 d 3 3 arcsec. This is in good agreement with the the ~ 12 arcsec 
radial extent of the core X-ray emission in the NS direction, suggesting that the compact nebula 
surrounding the pulsar is bounded by the pulsar wind termination shock. It is possible that 
this is actually a toroidal structure, much like that seen in the Crab Nebula, and that the elon- 
gated surface brightness distribution is the result of the inclination angle. In this interpretation, 
the axis of the toroid, which presumably lies along the rotation axis of the pulsar, lies in the 
east-west direction when projected onto the sky. We note that the long axis of 3C58 itself, as 
well as an extended jet-like feature shown in Figure 5 (right), are both aligned in this direction. 
Assuming that the radio filament lies along one side of the torus, its separation from the pulsar 




Figure 6: Left: ACIS image of 3C 58. Some portions of the nebula extend beyond the detector boundary, as 
indicated by the outermost contour from the HRC image, which is superposed. Center: ACIS image of the central 
region of 3C 58, with contours from 20 cm VLA data showing a faint radio wisp that bounds the X-ray core. 
Right: Saturated image of 3C 58 core, on a larger scale, revealing the faint jetlike feature extending toward the 
west. The dashed ellipse indicates the rough outline of the extended X-ray core. 



(~ 4.5 arcsec) implies an inclination angle of ~ 70°. 



4 Summary 

The strong shocks in SNRs and pulsar wind nebulae have long been regarded as sites where 
particles are accelerated to extremely high energies. SNRs, in particular, are leading candidates 
as the source of cosmic rays up to the knee of the spectrum. X-ray emission from PWNe 
also require particles with higher energies than expected in the postshock flow of the pulsar 
wind, indicating efficient acceleration at the termination shock. X-ray observations have now 
begun to provide direct evidence of the energetic particles and shock structures where this 
such acceleration takes place. Through studies of the dynamics and nonthermal X-ray emission 
from the shells of SNRs, and of the location and structure of the wind termination shocks and 
associated particle outflows in PWNe, strong constraints are being placed on models of the 
particle acceleration process. High resolution X-ray observations promise continued advances 
in this area, through measurements of SNR nonthermal emission and expansion rates, and the 
inner structure of PWNe. In addition, 7-ray observations with current and upcoming Cerenkov 
telescopes, as well as future space-borne observatories, hold considerable promise for probing 
these sites of extremely energetic particles. 
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